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Abstract 

Since the first experimental results were published in the 1990s, it has been believed that 
the sonoluminescence flash always occurs no more than a few nanoseconds before the minimum 
radius of a collapsing bubble. A concurrent belief has been that the period between sonolumi- 
nescence flashes is steady on the order of a few nanoseconds, and that sonoluminescence flashes 
occur with a "clock-like" regularity. To the contrary, data presented here show that the sono- 
luminescence flash can occur hundreds of nanoseconds before the minimum radius and that the 
sonoluminescence flash-to-flash period can vary on the order of hundreds of nanoseconds. These 
new findings may require a reexamination of the physics of sonoluminescence. 

1 Introduction 

Sonoluminescence is the emission of light flashes by sonically oscillating gas bubbles in a fluid. Sono- 
luminescence is brighter if a noble gas has been dissolved into the fluid - typically a water-glycerine 
mixture which has been degassed and then regassed with either argon, krypton or xenon. 

Our work was limited to an examination of single-bubble sonoluminescence, SBSL, because it is 
stable, and hence easier to observe and measure than multi-bubble sonoluminescence, MBSL - although 



we expect that what we have learned applies to MBSL as well. 

In SBSL, the sonoluminescent flash occurs only once per acoustic period, is very brief, and is very 
close to the moment of minimum bubble radius in the cycle of bubble expansion and contraction. 
Exactly how close to the time of minimum radius, and the degree to which that timing can vary, is 
the object of this investigation. We were motivated to repeat the timing studies done by others pQ, 
[5], not because we believed their results to be inconclusive, or unclear, but because we suspected 
that the observation of pulses with highly-regular timing might be an exceptional case resulting from 
a very well-controlled, smooth resonance, or some other conditions. We suspected that it might be 
possible to observe sonoluminescence in which the timing of the pulses was variable or farther away 
from the minimum radius. In fact, other published studies have hinted at this [3]. 

We have discovered that for the 23.55 kHz system which we constructed, the SL flashes can 
occur hundreds of nanoseconds before the time of minimum bubble radius. At those times, the fluid- 
dynamic models imply that the contents of the bubble should be close to room temperature, and that 
the motion of the bubble wall should be subsonic - far too slow to have created a shock-wave plasma, 
and well before any compressive heating could have occured. 
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2 Experimental Design 

To observe the motion of the bubble and the timing of the flash in the bubble cycle, we designed a 



Mie-scattering experiment (Figure 1 on page 14), similar to that used by other authors, [T], [3]. 

The procedure begins by preparing a mixture of water and 17% glycerine by volume, which has 
been degased and then brought into equilibrium with xenon, krypton or argon at 1/3 rd atmosphere. 
This solution is then poured into a 200 mL spherical glass flask with a narrow neck which has a 
resonant frequency around 23.55 kHz. 

A frequency generator, audio amplifier and impedance-matching inductor are used to drive PZTs 
attached to the sides of the flask, which is mounted on a labstand. When the frequency and drive are 
appropriate, bubbles can be trapped and made to glow by disturbing the surface of the water with a 
syringe. 

To begin the light scattering measurements, a 532 nm continuous-wave laser is aimed at the 
glowing bubble. The SL light and scattered laser light from the bubble are then focused onto a 
shielded Hamamatsu PMT, model R329-01, with a 3-inch condenser lens. The signal from the scope 
was recorded with a Tektronix 7104 digital sampling oscilloscope with a bandwidth of 3 GHz and a 
maximum sampling rate of lOGS/s. Our 2008 data were recorded with a 300 MHz Tektronix TDS 
3034 sampling oscilloscope with a maximum sampling rate of 3GS/s. This setup can be used to 
observe only the SL flash by turning off the laser. 
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3 Results 



3.1 Flash-to-Flash Jitter 

Before we ask how the timing of the flash varies with respect to the minimum bubble radius, let us 
first investigate how the timing between SL flashes varies. We call this the flash-to-flash jitter. 

We took these data with the Tektronix 7104 oscilloscope in November 2009. We recorded twenty 
scope traces of five sequential SL flashes like those shown in Figures [2] and |4j Thus we have a dataset 
of 100 flashes, from which we are be able to calculate 80 datapoints of flash-to-flash jitter. These were 
all taken with a resolution of 800 ps per sample. The flash times are identified with a simple algorithm 
which picks out the sample time corresponding to the maximum of the peak, (most negative voltage). 

Tables [T] and [2] summarize the flash data from Figures [2] and |4| We have presented the average 
period for the five flashes, the average frequency, the time of each of the five flashes, and the jitter, 
which is the amount which the four sequential flash-to-flash times differ from the average flash-to-flash 
time. The flash-to-flash jitter data from all 20 traces was calculated in this way and is presented as 
a histogram in Figure [5] Figure [3] shows a zoom-in on the first flash from Figure [2] (Note, due to 
constraints of space, we were not able to provide zoom-ins on each of the flashes from Figur^2| we 
will include these at a future time if we receive editorial permission.) 

Figure [2] and Table [T] show an exceptional case where the flash-to-flash jitter was as large as 102 

ns. 

In Figures [2] and [4] an oscillatory noise pattern with the same period as the SL is visible on the 
scope traces. It is more prominent in Figure [2] This noise is RF interference from the audio amplifier 
and cabling which induced a small voltage in the PMT signal. This noise signal has a gentle slope 
compared to sharpness of the SL flash, and should not affect our measurement of the peak SL time. In 
addition, the data in Figure [2] and [4] have a noise floor. (It is more prominent in Figure [2]). However, 
the SL signal is discernable above this noise floor, and because of the narrowness of the SL flash, this 
noise cannot move the peak SL time by more than 5 ns. 

Now, looking at Figure [3j we see that the SL flash has a full-width half-maximum of 5 ns, so that 
at most 10 ns of the jitter might be due to noise or the irregular response of the PMT to the flash. 
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However, the observations of 20 to 100 ns jitter as shown in the histogram in Figure [5] must be real 
and not due to systematic error. It has been previously reported that the SL flash has a "clock-like" 
precision, with a jitter on the order of picoseconds [BJ. 

Also note from examining Figures [2] and [4] how the intensity each flash can vary. 
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3.2 Flash-to-Minimum-Radius Timing 

Here we present our measurements of the timing of the sonoluminescent flash with respect to the 
minimum bubble radius. We performed this experiment in January and February 2008 and then 
again in November 2009. We present here a just a few of the hundreds of recorded flash events from 
sonoluminescing argon, krypton and xenon bubbles. We believe that to provide the most concise 
presentation and to eliminate the possibility of a phase offset from different signal paths, the laser 
scattering and SL flash should be observed with a single PMT. 

Figures [6] to 12 depict scope traces showing both the laser scattering and SL flash, and give an 



indication of our measurement of the time difference between the occurrence of the flash and the 
time of minimum bubble radius. This time difference can clearly be seen by eye, but to improve the 
likelihood that our estimate of the time difference is accurate we have also fit these scope traces to 
solutions of the Rayleigh-Plesset-Keller equation. 

Let us present a brief discussion of this equation and how its solutions can be fit to our raw laser 
scattering data. 

The oscillation of a gas bubble in a fluid can be described by the Rayleigh-Plesset-Keller equation: 



.pRR+ 3 -pR 2 (l-^- 



Here, R(t) is the bubble radius as a function of time, v s the speed of sound in the fluid, 77 the dynamic 
viscosity of the fluid, a the surface tension, P the ambient pressure, P(t) the pressure of the acoustic 
drive, and P gas the pressure of the gas in the bubble. 

This differential equation must be supplemented with an expression for P gas as a function of bubble 
radius: 

PUR) = PUV(R)) = { y^ h)1 (2) 

Here, b is the van der Waals excluded volume per mole. Such a formulation also requires that the 
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temperature of the gas in the bubble varies as: 



TqV 



(3) 



(V - nbY'- 1 



Note that a radial expansion or compression factor of ten will result in a hundred-fold decrease or 
increase in temperature, 7 being 5/3 for a monatomic gas such as xenon. 

The solution to the RPK equation Q is characterized by two experimental parameters - the 
equilibrium radius of the bubble Ro and the amplitude of the drive, Ro is the radius which the 
bubble would have if there were no drive present, which we can consider as the size of the injected 
bubble. When the bubble passes through or returns to this radius during its adiabatic cycle, it is at 
atmospheric pressure and room temperature. Since the RPK equation is second order in time, its 
solution also requires the specification of the initial velocity of the bubble wall, R(t = 0). However, 
these various solutions all settle into the same phase and response with respect to the drive after one 
or two periods have elapsed. So once the two parameters Ro and Pd are specified, a unique solution 
to the RPK equation, up to an overall time offset, can be generated. 

To match solutions of the RPK equation to light scattering curves requires an additional hypothesis 
on how the bubble scatters light and how the PMT responds to that light. Other researchers [T] have 
had success with the the hypothesis that the intensity of scattered light, and hence PMT voltage, is 
proportional to the cross sectional area of the bubble. Thus, 



where Vbkg is the amplitude of the laser light (PMT voltage) when no bubble is present. In order 
to calibrate the data we must find a value of a and Vbkg- Our primary goal in fitting the data is to 
obtain a value of t ffset, so that the time of minimum bubble radius can be determined from the fit. 

With this simple light scattering model, five parameters are required to fit voltage data to the 
solution R{t) of the RPK equation: R , P d , the time offset t ffset, the factor a and the background 
voltage Vbkg- Hence there will be a five parameter space in which to minimize the sum of squared 
differences. This can be challenging and unreliable, as the space may be filled with spurious local 




(4) 
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minima. Informed initial guesses of the parameters to be varied will improve the likelihood of good 
convergence. 

In practice, for our 23.55 kHz system, we have found that the values of Rq and Pd do not vary 
greatly and are typically around 9 fim and 1.5 aim respectively. Knowing this, we can fix i?o and Pd, 
and make the fitting procedure less computationally intensive. We performed all of these fits using 
the NonlinearModelFit routine in Mathematica v. 7.0. 

The data in Figure [6] depict two entire oscillation cycles. We can see the bubble grow, and then 
bounce several times before the cycle repeats. This data is fit to a solution of the RPK equation with 
i?o = 8.7 fim and Pd = 1.62 aim. At this scale the SL flashes are barely visible; but the next two 
Figures, [7] and |8j show zoom-ins on the first and second flashes, with the flash clearly visible before 
the minimum radius. In these two figures the fit was recalculated only allowing a, t offset, and Vbk g to 
vary to improve alignment with the minimum radius using a 2000 ns window of data. 

Figures [9] thru 12 depict zoom-ins near the minimum radius from the 2008 data set. These have 
been fit to solutions of the RPK equation with with Rq fixed at 9.5 fim, Pd fixed at 1.4 atm. Figure 



12 show a case where the peak of the flash is apparently 60 ns after the minimum radius. A close 



inspection of this figure suggests that the flash may have begun right at, or just a few nanoseconds 
after the minimum radius. If this is not truly a late flash, then at least it underscores the variability 
in flash timing when compared to the other traces. 

An added bonus of performing these fits is that we are able to obtain a calibration of the laser 
scattering data, and hence determine the velocity of the bubble wall at the moment of SL. Table [3] 
summarizes the timing and bubble wall velocity measurements from Figures [6] to 12 
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4 Conclusion 



The assertions that sonoluminescence has a "clock-like" regularity, or that the flash always occurs 
within a few nanoseconds of the minimum radius can be disproven simply by observing a single 
counterexample. We have presented several direct counterexamples to those assertions. 

We have attempted to depict our data in a manner which makes it obvious to the eye that the 
sonoluminescence flash timing is variable. It may be possible that a careful reviewer may find fault or 
take exception with some aspect of our method of fitting the Rayleigh-Plesset-Keller solutions to the 
light scattering data. Ideally, in a longer article, we would have room to go into more detail about 
our fitting methods. But we believe that the fits at least provide a corroboration of what is already 
quite apparent to the eye. If our calibration of the curves is accurate, then we have also shown the 
that bubble wall motion is subsonic at the moment of SL. 

It is obvious from our data that the sonoluminescence flash timing can vary on a scale much larger 
than previously reported. The sonoluminescent flash can occur hundreds of nanoseconds before the 
minimum radius, at times for which the bubble wall motion is subsonic, so that it seems that the flash 
cannot be the result of a thermal shock plasma formed by the collapse of the bubble. Thermal-shock 
models logically require that the flash occur right at or slightly before the minimum radius. For 
example, a shockwave model proposed by Greenspan [12J and Wu [8] requires the SL flash to occur 
just one tenth of a nanosecond before the minimum radius, because it is not until that moment in 
time that the RPK model implies the bubble wall velocity is large enough to create a thermal shock. 
Thus, that model and others like it must be ruled out. Furthermore, those models assume that the 
RPK model is an accurate description of bubble motion at those times, which may not be the case - 
those models do not take into account the boiling of the bubble wall that might dissipate the energy 
of the collapse before a thermal shock could be created. 

In a follow-up to this article we will introduce a new model of sonoluminescence which proposes 
that the SL flash is emitted by an excited cold condensate formed during the adiabatic expansive 
cooling which breaks apart to emit light and energetic electrons at random times near the minimum 
radius. 
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Table 1: A table of the five flash times from Figure [2j 



Flash 


Avg. Period (ns) Frequency (Hz) 


Flash Time ( 


ns) 


Jitter (ns) 


A 




8970.4 






B 




51405.6 




-27. 


C 


42462.2 23550.4 


93946.4 




78.6 


D 




136306. 




-102.2 


E 




178819. 




50.6 




Table 2: A table of the five flash times from Fij 


^ure 


l : 1 1 

4| 


Flash 


Avg. Period (ns) Frequency (Hz) 


Flash Time ( 


ns) 


Jitter (ns) 


A 




14980.8 






B 




57476.8 




-13.2 


C 


42509.2 23524.3 


99991.2 




5.2 


D 




142498. 




-2.0 


E 




185018. 




10.0 
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Table 3: A table of the times of the flash with respect to the minimum radius, At and bubble 
wall velocities from Figures [7] thru 12 Positive At values indicate time before the minimum radius, 
negative values, time after the minimum radius. 



Flash Figure At Bubble Wall Velocity 



Fi 


gure 


7 


63 ns 


-114 m/s 


Fi 


gure 


8 


47 ns 


-136 m/s 


Fi 


^ure 


9 


92 ns 


-91 m/s 


Fig 


ure 


10 


170 ns 


-62 m/s 


Fig 


ure 


11 


229 ns 


-51 m/s 


Fig 


ure 


12 


- 60 ns 


47 m/s 
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Figure Captions 

Figure [TJ An overhead schematic of the Flash Timing Experiment, showing approximate angles and 
distances. 

Figure [2j A scope trace showing five sequential SL flashes from a sonoluminescing krypton bubble. 
Figure [3} A zoom-in on the first flash from Figure [2] The peak of this flash is at 8970.4 ns. 
Figure [4j Another scope trace showing five sequential SL flashes from a sonoluminescing krypton 
bubble. 

Figure [5} A histogram of 80 flash-to-flash jitter measurements made from 100 krypton flashes. 
Figure [6[ A scope trace of a sonoluminescing krypton bubble with fit to the Rayleigh-Plesset-Keller 
model. This is from the 2009 dataset. 

Figure [7} A zoom-in near the first minimum radius of the data from Figure |6j showing the flash 
occuring about 63 ns before the minimum radius. The minimum radius time as determined from the 
fit is in approximate agreement with that determined by an "eyeball" reckoning. 
Figure [8} A zoom-in near the second minimum radius of the data from Figure [6j showing a flash 
about 47 ns before the minimum radius. 

Figure [9} A trace of a sonoluminescent argon bubble along with a three parameter fit to the RPK 
equation, with R fixed at 9.5 fim, Pd fixed at 1.4 atm and a, V^ g and t ffset allowed to vary. The 
results of the fit indicate that the flash is approximately 92 ns before the minimum radius. 



Figure 10 , A trace of a sonoluminescent argon bubble along with a three parameter fit to the RPK 
equation, with R fixed at 9.5 /im, Pd fixed at 1.4 atm and a, V^kg and t ffset allowed to vary. The 
results of the fit indicate that the flash is approximately 170 ns before the minimum radius. 
Figure A trace of a sonoluminescent krypton bubble along with a three parameter fit to the 
RPK equation, with Rq fixed at 9.5 jum, Pd fixed at 1.4 atm and a, Vbkg and t ffset allowed to vary. 
The results of the fit indicate that the flash is approximately 229 ns before the minimum radius. 



Figure 12, A xenon SL flash from the 2008 dataset . These data show the peak of the flash occuring 



60 ns after the minimum radius. 



13 




Figure 1: An overhead schematic of the Flash Timing Experiment, showing approximate angles and 
distances. 
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Figure 3: A zoom-in on the first flash from Figure [2] The peak of this flash is at 8970.4 ns. 
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Figure 4: Another scope trace showing five sequential SL flashes from a sonoluminescing krypton 
bubble. 



17 



i 




18 



o 
o 

CM 

o 



o 



J 

■ 

■ i 



• I 
• •I 
■ .it 



1 1 " 



ii. 

mil*! 











*U*n M|( 



■ t 

■ t 




mtittttt^ 



ititi+u'^ 




i: 



!■ 



* 1*444444 



■ i 
■1 








I ■ 



- 

- 

i 



--■ 



- 

- 
I 



s 

- 

I 




x 



o 
o 



o 
o 



o 

X 

o 

1 

o 



rT3 



rO 

fl 
O 

>, 
!h 

bJO 

.5 
'o 

X 

o 



j3 
o 

o 

x 



o 

o 
o 

X 

< 



CO 

Sh X 



19 



2 s 



o 



£ X) 

O 




O 
Sh 
2 



20 




21 



Pi 



t3 en 



o 
o 

< 



p cm 



I 

1 

o 



22 



■ ■ 




■V 



|£ 

<^ 03 
1* ^ 

co 
cu i — i 

cm ^ 
o « 

.-a x 

fn 

H rv 
O 

la 
*> co 

PL 05 

2 3 

CD -t-s 

-i-= 

o 43 

+J cu 



CU >iH 

I =s 

CO rj 

Qh +3 

CU Hh 

CU O 

9 CD 

* Go 



*S cu 

Pi 



o 
o 



O 

CU 

fl £ 

o o 
50 ^ 

fn co 

Pi (H 
CU o 

o 

CD _i 

s 1 

H h4 

O > 
Pi 

CD ° 
cd 

^ CO CD 

CO rrt CO 
fH M fH 

^ i i 

o ^ Pi 

1-1 V, " s 

cu cu 

fct a s 

bo cu 



23 



" 




•.; • > ' : 



cp 

x B 



8 

9 O 



pT ^ 

H CO 
.2 c6 

cm 

CO 

j5 cp 
cp 

W la 

^ CP 

CP 

£ cS 
Pi 



CP CP 

-H .-Pi 

CP 73 



O 

pi 

8 cS 
+3 cp 

CO r. 



bC > 

§ O 

cp cu 

3 I 



rH CD 

o m 

be 



Pi 

o 



.a a 



Pi ^ 
+^ Pi 

^ a 



X 



cp 

<n cm 



3 a 

h cs o) 



o 
o 
o 

CO 
Pi 
05 
CN 



24 




25 



